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Abstract 



In this paper we study dynamical chiral symmetry breaking in dimen- 
sionally regularized quenched QED within the context of Dyson-Schwinger 
equations. In D < 4 dimensions the theory has solutions which exhibit chiral 
symmetry breaking for all values of the coupling. To begin with, we study this 
phenomenon both numerically and, with some approximations, analytically 
within the rainbow approximation in the Landau gauge. In particular, we 
discuss how to extract the critical coupling Oc = ^ relevant in 4 dimensions 
from the D dimensional theory. We further present analytic results for the 
chirally symmetric solution obtained with the Curtis-Pennington vertex as 
well as numerical results for solutions exhibiting chiral symmetry breaking. 
For these we demonstrate that, using dimensional regularization, the extrac- 
tion of the critical coupling relevant for this vertex is feasible. Initial results 
for this critical coupling are in agreement with cut-off based work within the 
currently achievable numerical precision. 



I. INTRODUCTION 



It is fairly well established that quantum electrodynamics (QED), and in particular 
quenched QED, breaks chiral symmetry for sufficiently large couplings. This phenomenon 
has been observed both in lattice simulations as well as various studies based on the use 
of Dyson-Schwinger equations These latter calculations have generally relied on the 
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use of a cut-off in euclidean momentum in order to regulate divergent integrals, a procedure 
which breaks the gauge invariance of the theory. 

On the other hand, continuation of gauge theories to D < 4 dimensions has long been 
used as an efficient way to regularize perturbation theory without violating gauge invari- 
ance. In nonperturbative calculations, however, the use of this method of regularization is 
rarely used 0. Within the context of the Dyson-Schwinger equations (DSEs) only a few 
publications [^,0 have employed dimensional regularization instead of the usual momentum 
cut-off. 

It is the purpose of the present paper to study dynamical chiral symmetry breaking and 
the chiral limit within dimensionally regularized quenched QED. We are motivated to do 
this by the wish to avoid some gauge ambiguities occurring in cut-off based work, which we 
discuss in Sec. |T[ In that Section we also outline some general results which one expects 
to be valid for D < 4, independently of the particular vertex which one uses as an input 
to the DSEs. Having done this we proceed, in Section |T|, with a study of chiral symmetry 
breaking in the popular, but gauge non-covariant, rainbow approximation. Just as in cut-off 
regularized work, the rainbow approximation provides a very good qualitative guide to what 
to expect for more realistic vertices and has the considerable advantage that, with certain 
additional approximations, one may obtain analytical results. We check numerically that 
the additional approximations made are in fact quite justified. Indeed, it is very fortunate 
that it is possible to obtain this analytic insight into the pattern of chiral symmetry breaking 
in D dimensions as it provides us with a well defined procedure for extracting the critical 
coupling of the 4 dimensional theory with more complicated vertices. We proceed to the 
Curtis-Pennington (CP) vertex in Section 0. There we derive, for solutions which do not 
break chiral symmetry, an integral representation for the exact wavefunction renormalization 
function Z in D dimensions. We also provide an approximate, but explicit, expression for 
this quantity. The latter is quite useful, in the ultraviolet region, even if dynamical chiral 
symmetry breaking takes place as it provides a welcome check for the numerical investigation 
of chiral symmetry breaking with the CP vertex with which we conclude that section. Finally, 
in Section 0, we summarize our results and conclude. 



II. MOTIVATION AND GENERAL CONSIDERATIONS 

Although chiral symmetry breaking appears to be universally observed independently 
of the precise nature of the vertex used in DSE studies, it has also been recognized for a 
long time that the critical couplings with almost allQ of these vertices show a gauge depen- 
dence which should not be present for a physical quantity. With a bare vertex this is not 
surprising as this vertex Ansatz breaks the Ward-Takahashi identity. However, even with 
the Curtis-Pennington (CP) vertex, which does not violate this identity and additionally is 
constrained by the requirement of perturbative multiplicative renormalizability, a residual 



^Some vertex Ansatze exist which lead to critical couplings which are strictly gauge indepen- 
dent |^-1C|. However, these involve either vertices which have unphysical singularities or ensure 
gauge independence of the critical coupling by explicit construction. 
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gauge dependence remains [11,12 . 



Apart from possible deficiencies of the vertex, which we do not investigate in this paper, 
the use of cut-off regular ization explicitly breaks the gauge symmetry even as the cut- 
off is taken to infinity. This is well known in perturbation theory (see, for example, the 
discussion of the axial anomaly in Sect. 19.2 of Ref. [0) and was pointed out by Roberts 
and collaborators [|l^] in the present context. The latter authors proposed a prescription 
for dealing with this ambiguity which ensures that the regularization does not violate the 
Ward-Takahashi identity. 

As may be observed in Fig. |I], this ambiguity has a strong effect on the value of the 
critical coupling of the theory. The two curves in that figure correspond to the critical 
coupling ac of Ref. [l^ as well as the coupling one obtains by following the prescription 



of Roberts et al. It is straightforward to show, following the analysis of Ref. |]T2|], that these 
couplings are related through 

ol^ = . (1) 



Also plotted in this figure are previously published numerical results ||ll| , p!5i obtained with 



both of the above prescriptions. Note that, curiously, the critical couplings obtained with 



the prescription of Ref. fl^ (i.e. the calculation which restores the Ward-Takahashi identity) 
exhibits a stronger gauge dependence, at least for the range of gauge parameters shown in 
Fig. |. 

Gauge ambiguities such as the one outlined above are absent if one does not break the 
gauge invariance of the theory through the regularization procedure. Hence, we now turn to 
dimensionally regularized (quenched) QED. The Minkowski space fermion propagator S'(p) 
is defined in the usual way through the dimensionless wavefunction renormalization function 
ZijP') and the dimensionful mass function M{p^)^ i.e. 

The dependence of Z and M on the dimensionality of the space is not explicitly indicated 
here. Furthermore, note that to a large extent we shall be dealing only with the regu- 
larized theory without imposing a renormalization procedure, as renormalization |ll5| , p!6[| is 
inessential to our discussion. 

In addition to the above, we shall consider the theory without explicit chiral symmetry 
breaking (i.e. zero bare mass). This theory would not contain a mass scale were it not for 
the usual arbitrary scale (which we denote by v) introduced in D = 4 — 2e dimensions which 
provides the connection between the dimensionful coupling ao and the usual dimensionless 
coupling constant a = e^/47r: 

«£) = a i/^^ . (3) 

As u is the only mass scale in the problem, and as the coupling always appears in the above 
combination with this scale, on dimensional grounds alone the mass function must be of the 
form 

M(p2) = iya^^M(^,e] (4) 
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where M is a dimensionless function and in particular 

M(0) = i/a^ M(0,e) . (5) 

Moreover, as e goes to zero the v dependence on the right hand side must disappear and 
hence the dynamical mass M(0) is either zero (i.e. no symmetry breaking) or goes to infinity 
in this limit. This situation is analogous to what happens in cut-off regularized theory, where 
the scale parameter is the cut-off itself and the mass is proportional it. 

Note that M(0, e) is not dependent on a. This implies immediately that there can be 
no non-zero critical coupling in D 7^ 4 dimensions: if M(0) is non-zero for some coupling a 
then it must be non-zero for all couplings. 

Given these general considerations (which are of course independent of the particular 
Ansatz for the vertex) it behooves one to ask how this situation can be reconciled with a 
critical coupling ac of order 1 in four dimensions. In order to see how this might arise, we 
shall extract a convenient numerical factor out of M and suggestively re-write the dynamical 
mass as 

j_ 

M(0) = I. f— M(0,e) . (6) 

At present there is no difference in content between Eq. (^ and Eq. (|^). However, if we now 

define ac by demanding that the behaviour of M(0) is dominated by the factor ( as 
e goes to zero, which is equivalent to demanding that 

[M(0,6)]^ 1 , (7) 

then the intent becomes clear: even though M(0) may be nonzero for all couplings in D < 4 
dimensions, in the limit that e goes to zero we obtain 



M(0) j-^ a < ac (8) 

M(0) — > 00 a > ac . 



Note that in the above we have not addressed the issue of whether or not there actually 
is a M (0, e) with the property of Eq. (0). In fact, the numerical and analytical work in 
the following sections is largely concerned with finding this function and hence determining 
whether or not chiral symmetry is indeed broken for D < 4.| Notwithstanding this, as 
one knows from cut-off based work that there actually is a non-zero critical coupling for 
-D = 4, one can at this stage already come to the conclusion that M (0, e) exists and hence 
that quenched QED in D < 4 dimensions has a chiral symmetry breaking solution for all 
couplings. 

In summary, as the trivial solution M(p^) = always exists as well, we see that in 
D < 4 dimensions the trivial and symmetry breaking solutions bifurcate at a = while 



^The reader will note that as neither M(0, e) nor M (0, e) are functions of the coupling a, the 
value of Oc can be determined independently of the strength a of the self-interactions in D < 4 
dimensions. 
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for D = A the point of bifurcation is at a = a^, i.e., there is a discontinuous change in 
the point of bifurcation. As D approaches four (i.e. as e approaches 0) the generated mass 
M(0) decreases (grows) roughly hke (a/oc)^''^'' for a ^ Oc (~ Oc), respectively, becoming an 
infinite step function at a = Oc when e goes to zero. 



III. THE RAINBOW APPROXIMATION 



Let us now consider an explicit vertex. To begin with, we consider the rainbow approx- 
imation to the Euclidean mass function of quenched QED with zero bare mass in Landau 
gauge. It is given by 



M(p2 



(ei' 



e\2 



(3 - 2e) / 



1 



(27r)^ A;2 + M2(A;2) {p ~ kf 



(9) 



Note that the Dirac part of the self-energy is equal to zero in the Landau gauge in rainbow 
approximation even in D < 4 dimensions and hence that Z{p'^) = 1 for all p^. 

It is of course possible to find the solution to Eq. (|^) numerically - indeed we shall do 
so - however it is far more instructive to first try to make some reasonable approximations 
in order to be able to analyze it analytically. First, as the angular integrals involved in 
D-dimensional integration are standard (see, for example, Refs. and [jl^) we may reduce 
Eq. (^ to a one- dimensional integral, namely 



dk\ey-'M{e) 

p + M2(P) 



{l,e;2-e; ^] e{p' ~ k') 



+ -F\l,e-2-e-^^\e{k^~p^ 



where 



3 - 2e 



(4 



TT 



r(2 



(10) 



(11) 



Note that for D = A the mass function in Eq. (|10|) reduces to the standard one in QED4. 

In D ^ A dimensions the hypergeometric functions in Eq. ( [T0| ) preclude a solution in 
closed form. However, note that these hypergeometric functions have a power expansion in 
e so that for small e one is not likely to go too far wrong by just replacing these by their 
e = (i.e. D = 4) limit. After all, the reason for choosing dimensional regularization in 
the first place is in order to regulate the integral, and this is achieved by the factor of A;~2^^ 
not the hypergeometric functions. In addition, this approximation also corresponds to just 
replacing the hypergeometric functions by their IR and UV limits, so that one might expect 
that even for larger e that the approximation is not too bad in these regions 0. 



^It is however possible to show that a hnearized version of Eq. (^) always has symmetry breaking 
solutions even without making this approximation of the angular integrals. We indicate how this 
may be done in Appendix A. 
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Making this replacement, i.e. 



k^ + M^{k^) 



k^ + M^{k^ 



(12) 



allows us to convert Eq. (JWj) into a differential equation, namely 



with the boundary conditions 



p2 + ]Vf2(p2 



-M{p') = 



p'M'{p%2^, = 0, 



p^M{p^ 



\p^=oo 



. 



(13) 



(14) 



Unfortunately, the differential equation (|T3|) still has no solutions in terms of known special 
functions. Since the mass function in the denominator of Eq. (12) serves primarily as an 
infrared regulator we shall make one last approximation and replace it by an infrared cut-off 
for the integral, which can be taken as a fixed value of M'^{k'^) in the infrared region (for 
convenience we shall call this value the 'dynamical mass' m). This simplifies the problem 
sufficiently to allow the derivation of an analytical solution. 

In terms of the dimensionless variables x = p^/z/^, y = k'^/v'^ and a = vr^jv^ the 
linearized equation becomes 



M{x) = ac^ 



dyy^-'M{y) , fdyy-'M{y) 



y 



+ 



y 



(15) 



for simplicity, we do not explicitly differentiate between M{x) and M{p^). This may be 
written in differential form as 



x^M\x)\ + ac, x-'M{x) = 
with the boundary conditions 

M'(x)U=, = 0, [xM(x)]' U=oo = 



. 



This differential equation has solutions in terms of Bessel functions 



M(x) = 



Ci J\ 



ex2 



+ C2 J-> 



ex2 



(16) 



(17) 



(18) 



where we have defined A = 1/e in order to avoid cumbersome indices on the Bessel functions. 
The ultraviolet boundary condition Eq. (p!?]) gives C2 = while the infrared boundary 
condition leads to 



+ 



ex2 



X2 



6X2 



. 



(19) 



This equation may be simplified using the relation among Bessel functions 
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zJ'^{z) + \J\{z) = zJx-i{z) 



(20) 



and becomes 



—Jx-i I — 



e X2 \ e X2 



. (21) 



Clearly this equation is satisfied by Ci = 0, which corresponds to the trivial chirally sym- 
metric solution M{x) = 0. However, for values of a which are such that the argument of 
the Bessel function in Eq. (|2T| ) corresponds to one of its zeroes, the equation is also satisfied 
for Ci 7^ 0, i.e. for these values of a there exist solutions with dynamically broken chiral 
symmetry. If we define Ja-i,! = \/ ^OLC^j ta"!"^ to be the smallest positive zero of Eq. (pTD, 
the dynamical mass for this solution becomes 

1 

m = V a^/'^ = va^- \ | . (22) 

\e jA-i,!/ 

Note that for this solution the normalization Ci is not fixed by Eq. ([ISD as this equation is 
linear in M{x). Later on we shall fix Ci by demanding that M(a) = m, however there is 
no compelling reason to do this and one might alternatively fix the normalization in such a 
way as to approximate the true (numerical) solutions of Eq. (^) as well as possible. Finally, 
note that, as expected, a dynamical symmetry breaking solution exists for any value of the 
coupling and that the expression for the dynamical mass is in agreement with the general 
form expected from dimensional considerations [i.e. Eq. (^]. 

In order to extract etc, we need to look at the behaviour of m as e goes to zero (i.e. 
A — >• oo). This may be done by noting that the positive roots of the Bessel function Jx have 
the following asymptotic behaviour (see, for example, Eq. 9.5.22 in Ref. |[l8|] ): 



JV-Az(C) + f:|S' C = A-^/^a„ (23) 

where is the sth negative zero of Airy function A\{z)^ and z{C) is determined (-2(C) > 1) 
from the equation 

-(-()3/2 = _ arccos -. (24) 

3 z 

For large A the variable C, is small and so it is valid to expand z around 1. Writing z = 1 + 5 
we obtain 

arccos ^ ~ - ^^3/2 (25) 
1 + 5 12 ^ ^ 

and so 5 ~ —C,/2^/^ yielding, to leading order, 

- = ■ (26) 

If we define 
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7 = ~ 1.855757 (27) 

then the leading terms in the expansion of Ja-i,i are 

jA-i,i~A + 7Ai/='-l + 0(A-i/3). (28) 
Also, the coefficient appearing in Eq. (E^ may be expanded 



c, + rf = ln(47r) + i + 7/;(l) (29) 

47]- J 



SO that for small e the dynamical mass becomes 

a + de) 



e^+i-i^-'" . (30) 



Note that the behaviour of the ffist term (for e going to zero) dominates over the exponential 
function, as required in Eq. (|^). Hence the critical coupling in four dimensions is given by 
7r/3, as expected from cut-off based work 

Returning now to the mass function itself, we may substitute the expression for the 
dynamical mass, i.e., Eq. (p2D, together with our choice of normalization condition 



M(p^ = m^) = m , (31) 
into Eq. (^) in order to eliminate Ci. One obtains 

\v\ J\\3\-iM 

Note that the explicit dependence on v (and hence a) has been completely replaced by m 
in this expression. 

So far we have taken a independent of the regularization. As we have seen this leads to 
a dynamically generated mass which becomes infinite as the regulator is removed. Fomin 
et al. [0 examined (within cut-off regularized QED) a different limit, namely one where the 
mass m is kept constant while the cut-off is removed. In our case this limit necessitates that 
the coupling a is dependent on e through 



a 



\{l + 2^ei) (33) 



[see Eq. (|30D; note that etc is approached from above]. The limit may be taken analytically 
in Eq. ( |5^ by making use of the known asymptotic behaviour of the Bessel functions (see 
Eq. 9.3.23 of Ref. [0), i.e. 

J,(\ + \'''z) ^ {^^"\\{-2'''z) , (34) 
as well as the asymptotic expansion of ja-i i in Eq- (11) • One obtains 
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M(p) = — (ln^ + 1) , (35) 

which agrees with the resuh in Ref. 0]. 

To conclude this section, we analyze the validity of the approximations made by solv- 
ing Eq. (|^) numerically and comparing it to the Bessel function solution in Eq. (p2D. In 
Fig. ^ we have plotted the mass function (divided by v) as a function of the dimensionless 
momentum x for a moderately large coupling [a = 0.6) and e = 0.03. The solid curve corre- 
sponds to the exact numerical result [Eq. @] while the dashed line is a plot of Eq. (|32|) for 
these parameters. As can be seen, the approximation is not too bad and could actually be 
made significantly better by adopting a different normalization condition to that in Eq. (|3TD 
However, no further insight is gained by doing this and we shall not pursue it further. 

One might naively think that most of the difference between the Bessel function and the 
exact numerical solution comes from the linearization of Eq. (^ - i.e., the approximation 
made by going from Eq. ([T^) to Eq. (|15]), as the only approximation made prior to this is 
to replace the hypergeometric functions by unity, which is expected to be good to order 
e (i.e. in this case, 3 %). This turns out to be not the case; the dotted curve in Fig. ^a 
corresponds to the (numerical) solution of Eq. ([T^) . Not only is the difference to the true 
solution essentially an order of magnitude larger than expected (about 30 % - note that 
Fig. ||a is a log- log plot), it is actually of opposite sign to the equivalent difference for the 
Bessel function. In other words, the validity of the two approximations is roughly of the 
same order of magnitude and they tend to compensate. 

Why are the quantitative differences rather larger than expected? On the level of the 
integrands the approximations are actually quite good. In Fig. we show the integrands 
of Eqs. (P), (p!5|) and Eq. (|12|) for a value of x in the infrared [x ~ 7.1 10^^^). Clearly the 
replacement of the hypergeometric functions by unity is indeed an excellent approximation, 
as is the linearization performed in Eq. (|l^) (except in the infrared, as expected. Note 
that when estimating the contribution to the integral from different y one should take into 
account that the x-axis in Fig. ^ is logarithmic). The real source of the 'relatively large' 
differences observed for the integrals in Fig. ^ is the fact that these are integral equations 
for the function M{x) - small differences in the integrands do not necessarily guarantee 
small differences in M{x). To illustrate this point, consider a hypothetical 'approximation' 
to Eq. (|T5|) in which we just scale the integrands by a constant factor 1 + e and ask the 
question how much this affects the solution M{x). For x = the answer is rather simple: 
the hypothetical approximation just corresponds to a rescaling of a by 1 + e and as M(0) 
scales like a 2^ we find that the solution has increased by a factor (1 + e) 27. In other words, 
even in the limit e — > there remains a remnant of the 'approximation', namely a rescaling 
of M(0) by a factor e^/^ ~ x gi 



IV. THE CURTIS-PENNINGTON VERTEX 

We shall now leave the rainbow approximation and turn to the CP vertex. The expres- 
sions for the scalar and Dirac self-energies for this vertex, using dimensional regularization 
and in an arbitrary gauge, have already been given in Ref. 0. Before we discuss chiral 
symmetry breaking for this vertex we shall first examine the chirally symmetric phase. We 
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remind the reader that in this phase in four dimensions the wavefunction renormahzation 
has a very simple form for this vertex [19|, namely 




Z{x,fi^)\M{x)=0 = 

where the renormalized Dirac propagator is given by 

Sip) = . 



(36) 



(37) 



Here ^ is the gauge parameter and yU^ is the (dimensionless) renormahzation scale. This 
power behaviour of Z{x) is in fact demanded by multiplicative renormalizability pO[ as well 
as gauge covariance We shall derive the form of this self-energy in Z) < 4 dimensions, 
which will provide a very useful check on the numerical results even if M{x) 7^ as long as 



X » 



' M(x) 



A. Z{p'^) in the chirally symmetric phase 

In the chirally symmetric phase, the unrenormalized Z[x) corresponding to the CP vertex 
in D dimensions is given by 



, — e 



47r r(2 — e) Jo x — y 



\xJJ v + x \x 



(38) 



This equation may be obtained from Eq. (A6) of Ref. 0] by setting b{y) equal to zero in 
that equation and by using Eq. (A8) of the same reference in order to eliminate the terms 
with coefficient a?{y)- The angular integral Ii{w) is defined to be 

/f (u;) = (1 + w) 2i^i(l,e;2 -e;w) Q<w<l (39) 
(«;) = If (t/7-i) w>l . (40) 



In four dimensions the solution to Eq. ( |38D is given by a Z[x) having a simple power 
behaviour while for D < A this is clearly no longer the case. Nevertheless, it is possible 
to derive an integral representation of the solution of Eq. ( pSj ) by making use of the gauge 
covariance of this equation. We do so in Appendix B, with the result 



Z(x) = 2^-^ r(2 - e) / c/MM^-V-™'V2-e(v^M) • (41) 

Jo 

Although this result is exact it is somewhat cumbersome to evaluate numerically because, 
for e — > 0, the oscillations in the integrand become increasingly important. For this reason 
we shall approximate the integrand in Eq. ( pH]) by its IR and UV limits, as we did for the 
rainbow approximation (as before, this approximation is good to order e). Using 

Jf (w;) = 1 + -^w + 0[w^] (42) 
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this approximation yields 

a {Any 



Z x) 



1 + 



47rr(2 



e Jo Jx 



(43) 



This may be converted to the differential equation 



Z"{x) 

where c is defined to be 



-Z'(x) 



X 



X 



1+e 



Z'(x) + 2 - — - Z(x) 



X 



^ ~ 2^r(3-e) 
and the appropriate boundary conditions are 

x^~^Z(x)\^.,=n = 0, Z(x 



1. 



(44) 



(45) 



(46) 



(The IR boundary condition arises from the requirement that the integral in Eq. (|43|) needs 
to converge at its lower limit.) In order to solve Eq. (0), it is convenient to change variables 
to 



X 



2-e 
and to define 

e 

so that the differential equation becomes 

zZ" — a{l — z)Z' — a{a — 1)Z 
while the boundary conditions now are 



z-^Zl 



0, 



z 



z=0 



1. 



(47) 



(48) 



(49) 



(50) 



This equation is essentially Kummer's Equation (see Eq. 13.1.1 of Ref. [|T8l; we use the 
notation of that reference in the following). Its general solution may be expressed in terms 
of confluent hypergeometric functions, i.e. 



Z = z^+^e-"' 



Ci M{a, a + 2; az) + C2U{a,a + 2; az) 
{Ci [7(0 + 1,-02;) + az'y{a,-az)] + C2[l + -2]} . 



(51) 



The UV boundary condition is fulfllled if C2 = 1 while Ci is not flxed by the boundary 
conditions. Although Eq. (|4^) is solved by Eq. (|5TD for arbitrary Ci we shall concentrate 
on the solution with Ci=0. The reason for this is that the solution to the unapproximated 
integral [Eq. (0)] vanishes at x = (see Appendix B) while the term multiplying Ci in 
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Eq. 
Eq. 



ID diverges like x^*^ ^ and is therefore unlikely to provide a good approximation to 
Hence we obtain 



Zix) 



1 + 



c 



-X 



( ^ - 

exp( — X 
e 



(52) 



2-e 

Finally, the renormalized function ^(x,/i^) may be obtained from this by demanding that 
Zijj?^ /i^) = 1 so that the renormalized wavefunction renormalization becomes 



1 + 



-x~ 



1 + 2^/^"^' exp( 



exp(— % ^) 
) 



(53) 



Only in the limit D — > 4 does this reduce to the usual power behaved function found in 
cut-off based work [Eq. (|36|)] while for D < 4 it vanishes non-analytically at x = 0. On the 
other hand, note that the solution to Eq. (|38|) - for finite e - only goes to zero linearly in x. 
For the purpose of this paper this difference in the analytic behaviour in the infrared does 
not concern us as for solutions which break chiral symmetry the infrared region is regulated 
by M^(x) so that we do not expect the chirally symmetric ^ to be a good approximation 
in this region in any case. 



B. Chiral symmetric breaking for the CP vertex 

We shall now examine dynamical chiral symmetry breaking for the CP vertex in the 
absence of any explicit symmetry breaking by a nonzero bare mass, as before. Even for 
solutions exhibiting dynamical symmetry breaking, it is to be expected that the analytic 
result derived for Z{x)\M{x)=a [Eq- O)] remains valid as long as x is large compared to 
{M{x)/vY and e is sufficiently small. That this is indeed the case is illustrated in Fig. |^, 
where we show a typical example of Z~^[x) for a solution which breaks chiral symmetry. In 
this Figure, as well as in the rest of this Section, we shall be dealing with the renormalized 
Z{x) and M{x) instead of the unrenormalized quantities in the previous sections. This 
makes no essential difference to the physics of chiral symmetry breaking, although it of 
course effects the absolute scale of Z{x). For a discussion of the renormalization of the 
dimensionally regularized theory we refer the reader to Ref. [0. 

The comparison to the analytic result in Fig. ^ provides a very convenient check on the 
numerics. Another check is provided by plotting the logarithm of M(0) against the logarithm 
of the coupling. According to Eq. (^ this should be a straight line with gradient ^. As can 
be seen in Fig. ^ not only does one observe chiral symmetry breaking down to couplings as 
small as a = 0.15, the expected linear behaviour is confirmed to quite high precision. 

Although the numerics in Z) < 4 dimensions are clearly under control, the extraction 
of the critical coupling (appropriate in four dimensions) has proven to be numerically quite 
difficult. From the discussion in Sections and |T|, we anticipate that the logarithm of the 
dynamical mass has the general form 
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where the last term is subleading as compared to the first as e tends to zero. For sufficiently 
small e, therefore, is related to the gradient of In (M(0)) plotted against e~^. 

In Fig. 1^ we attempt to extract in this way. The logarithm of M(0) was evaluated for 
e ranging from 0.03 down to e = 0.015 for a fixed gauge ^ = 0.25. The squares corresponds 
to a coupling constant a = 1.2, although some of the points at lower e have actually been 
calculated at smaller a and then rescaled according to Eq. (||). At present we are unable, for 
these parameters, to decrease e significantly further without a significant loss of numerical 
precision. (We also note in passing that it is quite difficult numerically to move away from 
small values of the gauge parameter; ^ = 20, which, judging from Fig. ^ would not require 
a very high numerical accuracy for ac, is unfortunately not an option.) 

The two fits shown in Fig. |^ correspond to two different assumptions for the functional 
form of M(0, e), which is a priori unknown. The curves do indeed appear to be well approxi- 
mated by a straight line, however we caution the reader that this does not allow an accurate 
determination of ac as the gradient is essentially determined by the 'trivial' dependence on 
log(Q;) (more on this below). The solid line corresponds to the assumption that the leading 
term in M(0, e) has the same form as what we found in the rainbow approximation, i.e. 

■■'(^)^^-(^)-'(^)* • (-) 

With this form the ffi parameters ac and ci are found to be 

ac = 0.966 ci = -1.15 . (56) 

Indeed, the critical coupling is similar to what is found in cut-off based work (see Sec. ||; 
in Ref. |TB| the value was 0.9208 for ^ = 0.25). At present it is difficult to make a more 



precise statement, let alone differentiate between the two curves plotted in Fig. |l|, as ac is 
quite strongly dependent on the functional form assumed in Eq. (|55|) . In fact, allowing an 
extra constant term on the right hand side of Eq. (|55|) reduces the critical coupling to 0.920 
and the addition of yet a further term proportional to es increases it again to 0.931. As these 
numbers appear to converge to something of the order of 0.92 or 0.93 one might think that 
ac has been determined to this precision. However, it is not clear that the functional form 
suggested by the rainbow approximation should be taken quite this seriously. The dashed 
line in Fig. § corresponds to a fit where the power of e of the subleading term has been left 
free, i.e. 

The optimum ffi assuming this form for In (M(0)) yields a power quite different to | and a 
very much smaller ac'- 

ac = 0.825 ci = -0.801 C2 = 0.688 . (58) 

To conclude this section, let us discuss why it is that the functional form of the subleading 
term M(0, e) appears to be rather important even if e is already rather small. The reason 
for this is two-fold: most importantly, although the leading e dependence of ln(M(0)) is 
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indeed e~^, the coefficient of tliis term (leaving out tfie trivial a dependence) is ln(Q;c)- 
As etc is rather close to 1 one therefore obtains a strong suppression of this leading term, 
increasing the relative importance of the subleading terms. In addition, it appears as if the 
numerical results favour a subleading term which is not as strongly suppressed (as a function 
of e) as suggested by the rainbow approximation (i.e. the power of of the subleading 
term appears to be closer to | rather than |). This again increases the importance of the 
subleading terms. 

V. CONCLUSIONS AND OUTLOOK 

The primary purpose of this paper was to explore the phenomenon of dynamical chi- 
ral symmetry breaking through the use of Dyson- Schwinger equations with a regularization 
scheme which does not break the gauge covariance of the theory, namely dimensional reg- 
ularization. It is necessary to do this as the cut-off based work leads to ambiguous results 
for the critical coupling of the theory precisely because of the lack of gauge covariance in 
those calculations. In particular, this should be kept in mind when using the expected gauge 
invariance of the critical coupling as a criterion for judging the suitability of a particular 
vertex. 

To begin with, we have shown on dimensional grounds alone and for an arbitrary vertex, 
that in < 4 dimensions either a symmetry breaking solution does not exist at all (in 
which case, however, it would also not exist in D = 4 dimensions) or it exists for all non 
zero values of the coupling (in which case a chiral symmetry breaking solution exists in D = 4 
for a > Qfc)- For Dyson- Schwinger analyses employing the rainbow and CP vertices we have 
shown that it is the second of these possibilities which is realized. For these symmetry 
breaking solutions the limit to D = 4 is necessarily discontinuous and so the extraction of 
the critical coupling of the theory (in 4 dimensions) is not as simple as in cut-off regularized 
work. 

We next turned to an examination of symmetry breaking in the rainbow approximation 
in Landau gauge, both analytically and numerically. Indeed, for this vertex one could rewrite 
the (linearized) Dyson-Schwinger equation as a Schrodinger equation in 4 dimensions and 
appeal to standard results from elementary quantum mechanics to explicitly show that the 
theory always breaks chiral symmetry if D < 4. We also showed how the usual critical 
coupling CKc = I may be extracted from the dimensionally regularized work. 

We concluded this work with an examination of the CP vertex. By making use of the 
gauge covariance of the theory we derived an exact integral expression for the wavefunction 
renormalization function Zijp') of the chirally symmetric solution. Furthermore we obtained 
a compact expression for this quantity which is an excellent approximation to the true 
Z{p^) even for solutions which break the chiral symmetry. Finally, we extracted the critical 
coupling corresponding to this vertex and found that, within errors, it agrees with the 
standard cut-off results. 

In the future, we plan to increase the numerical precision with which we can extract this 
critical couphng for the CP vertex by an order of magnitude or so. The factor limiting the 
precision at present is that when solving the propagator's Dyson-Schwinger equation with 
the CP vertex by iteration the rate of convergence decreases dramatically as e is decreased 
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below e ~ 0.015. If this increase in precision can be attained it will enable one to make a 
meaningful comparison with the cut-off based results shown in Fig. |l|. 
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APPENDIX A: CHIRAL SYMMETRY BREAKING IN RAINBOW 

APPROXIMATION 



In this appendix we show that the linearized version of Eq. (^), i.e. 

has symmetry breaking solutions for all values of the coupling. Our aim here is to convert 
this equation to a Schrodinger-like equation, which we do by introducing the function 

"^^'^ - J {2n)^ k^ + m^ • ^^^^ 

With this definition we have 

(-□ + m^) ^(r) = / ^e^''M{k') (A3) 
where □ is the D-dimensional Laplacian and so 

After shifting the integration variable {p ^ p + k) the last equation can be written in the 
form of a Schrodinger-like equation 

H%lj{r) = -m^V'(r), (A5) 

where H = —O -\- V (r) is the Hamiltonian, E = —m? plays the role of an energy and the 
potential V{r) given by 

where 
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Air 



2-e 



(A7) 



For D = 3 the coefficient i] is 2ua while near D = 4 it is ^az/^^ It is well known from any 
standard course of quantum mechanics (see, for example, Ref. ||2T|]) that potentials behaving 
as at infinity, with s < 2, always support bound states (actually, an infinite number of 
them). In the present case this can be seen by considering the Schrodinger equation ( |A5[ ) 
for zero energy, i.e. E = 0. The s-symmetric wave function then satisfies the equation 



r + 



The solution finite at the origin r = is 



(A8) 



(A9) 



The Bessel function in ( |A9| ) has an infinite number of zeros, which means that there is an 
infinite number of states with E < 0. 

Returning now to Eq. ( [A5D , we can estimate the lowest energy eigenvalue variationally 
by using 



ip(r) = const. X Ji^i { — i^'^ j ■ 



ipir) = Ce~ 



(AlO) 



as a trial wavefunction. Here C is related to k, by demanding that ijj is normalized, i.e 



\C\ 



(All) 



where O^) is the volume of a D-dimensional sphere. Calculating the expectation value of 
the "Hamiltonian" H on the trial wave function in Eq. ( |A10D we find 



E,(k^] 



\H\ 



)D-2 



T{D) 



(A12) 



The minimum of the "ground state energy" in Eq. ( |A12| ), Eo{k,), is reached at 



A-D 



{D-2) 



T{D) 



(for = 3 the parameter n is va while near D = A it is v 
expression 



1 



-m 



D 



tt/2 

,D-A 
D-2 



(A13) 

and is given by the 
(A14) 



where the 1 is the contribution from the kinetic energy while the — l) corresponds to 
the potential energy. For D > 2 the potential is attractive and for 2 < D < 4 it is always 
larger than the kinetic energy, so for this case we get dynamical symmetry breaking for any 
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value of a. For example, for D = 3, one obtains = —k^ = —v^o? which coincides precisely 
with the ground-state energy of the hydrogen atom (not surprisingly, as we have used the 
ground-state hydrogen wave function as our trial function). In this case the dynamical mass 
is m = va. 

For D near 4, on the other hand, we obtain from Eq. ([A14|) that 



a 



2e 



This is of the general form anticipated in Section with ac = f • Indeed, for D = 4, the 
Schrodinger equation (|A5|) becomes an equation with the singular potential 

y{r) = -\, r/ = 4-. (A16) 



r 



7r/3' 



Again, it is known from standard quantum mechanics that the spectrum of bound states 
for such a potential depends on the strength 77 of the potential: it has an infinite number 
of bound states with < if > 1 and bound states are absent if 77 < 1. Thus, the true 
critical value for the coupling is expected to be etc = ^r/S instead of the etc = 7r/2 obtained 
with the help of the variational method (which made use of the exponential Ansatz for the 
wavefunction and thus only gave an upper bound for the energy eigenvalue). 

APPENDIX B: CHIRALLY SYMMETRIC QED FROM THE 
LANDAU-KHALATNIKOV TRANSFORMATION 



Because the CP vertex in the chirally symmetric phase of QED is gauge-covariant it is 
possible to derive an integral representation of the wavefunction renormalization function 
Z{x) [see Eq. (|38|) ] from the Landau-Khalatnikov transformation |23|. This transformation 



relates the coordinate space propagator S^{u) in one gauge to the propagator in a different 
gauge. Specifically, with covariant gauge fixing, we have 

~S<^{u) = e4--'^I^(o)"^Wl5«=°(M) (Bl) 

where A(m) is essentially the Fourier transform of the gauge-dependent part of the photon 
propagator, i.e. 

A(u 



Specifically, we obtain 



where 



S^{u) = e-^^("")''5«=°(M) (B3) 



-^r(-6)(7r)^e ■ (B4) 



Substituting the coordinate-space propagator in Landau gauge, i.e. 
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u 



2 7rD/^ 



u 



D ' 



(B5) 



and carrying out the inverse Fourier transform of Eq. ( P3|) one obtains the wavefunction 
renormahzation function in an arbitrary gauge, namely 



Z{x) 



2 7r^/2 

x5 2^-^r(2 



D 



u 



(B6) 



Note that for small x this function vanishes: 

re) 



Z(x) 



4e (2 



r <= X + O ( a; 



(B7) 



It may be checked explicitly that Eq. ( p6|) is indeed a solution to Eq. (|38|) for arbitrary 
D by making use of the expansion of Eq. ( P6| ) around x"*^ = 0. To be more precise, consider 
the RHS of Eq. ( p8|) upon insertion of the power in the place of Z{y). Note that the 
integral converges only ife>(5>e — 2. After some work the result is that the R.H.S. of 
Eq. (pSD becomes 



1 + 5- 



2-e 



-X 



V{2-e) ^ 



r(2-e + 72)72-5 + 6 



(B8) 



For e < 1 this may simplified further by applying Dougall's formula (Eq. 1.4.1 in p4[) which, 
in this case, reduces to 



E 



r(e + n) 



1 



TT 



1 



r(2 - e + n) n - (5 + e sin(7re) sin(7r[e - 6]) V{1 - 5)T{2 + 5 - 2e) 
Using this result, Eq. (|B8|) becomes 



(B9) 



1 _ ,2 - 6 . , r(l - 6)r(2 - 6)r(6 - 5)V{2 + e - 5) 
2 r(2 + 5-2e)r(l-5) 



(BIO) 



Note that, as opposed to the integral representation Eq. (pS]), this expression is defined for 
5 outside the range e > 5 > e — 2 and so we may use it as an analytical continuation of the 
integral. Furthermore, note that this last expression vanishes for integer 5 > 1 hence we 
cannot obtain a simple power expansion around x = for Z{x) in this way. 

On the other hand, an expansion in powers of x~'^ is possible. If we seek a solution of 
the form 



z{x) = E 



Cri.X 



(Bll) 



n=0 
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we may equate the coefficients of equal powers of after inserting the series (|B11| ) into 
both sides of Eg . (|38|) . This way we obtain the recurrence relation for the coefficients c„ 
(co = 1) as 



Cn+l 



Cn 2(n+l, 
This may be solved leading to 



r(-6)r(3-e: 



r(l + en + e)r(2-e-en) 
r(2 - 2e - en)T{l + en) 



(B12) 



(B13) 



so that finally we obtain 
Z{x) = r(2 



rn + ne) 



ne)n\ 



gr(-e)r(3 - 6) 



-X 



(B14) 



as the series expansion of the solution to Eq. (|38D . The reader may check that this coin- 
cides precisely with the corresponding expansion of the solution obtained via the Landau- 
Khalatnikov transformations [Eq. (|B6| )]. The latter may be obtained by changing the vari- 
able of integration from u to u/ y^, expanding the exponential in the integrand and making 
use of the standard integral 



X 



2a-l^-a 



2 



r(i + 



v—a 
2 



(B15) 
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FIGURES 
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FIG. 1. The critical coupling for the CP vertex. The solid line is taken from the bifurcation 
analysis carried out in Ref. |12|, which agrees with the numerical results (open squares) of Ref. |11 



The dashed line corresponds to the bifurcation analysis carried out with the 'gauge violating term' 
removed (as suggested in Ref. [|14|) and agrees with the numerical results (open triangles) of 
Ref. [ITsi. 
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FIG. 2. The mass function for rainbow QED for a = 0.6 and e = 0.03 as a function of 
X = p^/u"^. The dynamical mass is specified by m/v = 2.24 x 10~^, where v is the scale intro- 
duced by dimensional regularization. The solid line corresponds to the exact numerical solution of 
Eq. (^) , the dashed line is the Bessel function of Eq. (^) and the dotted line is the solution of the 



-10 
log(y) 
0.6 and e 



Dyson- Schwinger equation with the hypergeometric function replaced by unity (Eq. 12). In (a) the 
actual mass function is shown, while in (b) we show the integrand at a particular value of x. 
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FIG. 3. A typical (inverse) wave function renormalization function Z ^{x,ii^) corresponding 
to a chiral symmetry breaking solution. Note that the mass function M{x) is of the same order as 



z- 



itself. Nevertheless, the analytical chirally symmetric solution of Sec. IV A| (thin line) 



provides an excellent approximation (better than one part in a thousand for x > fi ) for x > M{x). 
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FIG. 4. The logarithm of the dynamical mass as a function of ln(a) for e = 0.025 (i.e. ^ = 20). 
The gauge parameter is fixed at ^ = 0.25 and the renormalization point is = 10^. The open 
squares are the numerical values while the solid line is a linear fit to these points. Note that the 
dependence on the coupling expected in Eq. (|5|) is reproduced to high precision. 
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FIG. 5. The logarithm of the dynamical mass as a function of ^ for a coupling of a = 1.20. 
All other parameters are as in Fig. ^. The open squares are the numerical values while the two 
lines are fits (see main text) 
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